APPROXIMATE  SELF-CONSISTENT  FIELD  CALCULA< 
TfOMS  OF  MOLECULAR  POTENTIAL.  SURFACES 

CHANG  SOP  YUN 


will 

"  ■'ASM  \\M 

Mil 


■t  ! 


iflHHHi 


■:■  ^"rrr': }:.'.; 


APPROXIMATE 
SELF-CONSISTENT  FIELD  CALCULATIONS 
OF  MOLECULAR  POTENTIAL  SURFACES 

by 


YUN,  CHANG  SOP 

Lieutenant,  Republic  of  Korea  Navy 

B.S.,  Naval  Academy,  1961 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  CHEMISTRY 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
May  1967 


ABSTRACT 


The  Approximate  Self -Consistent  Field  Molecular  Orbital  Theory 
(ASCFT)  with  complete  neglect  of  differential  overlap  (CNDO)  has  been 
applied  to  the  calculation  of  electronic  structure  for  several  selected 
molecules  in  this  text.  The  theory  leads  to  calculated  equilibrium 
geometry  configurations  and  the  barrier  to  internal  rotation  of  ethane, 
inversion  barrier  of  ammonia  and  the  barrier  for  cis  to  trans  conversion 
for  difluorodiazine.  These  results  are  in  reasonable  agreement  with 
experimental  values  in  most  cases.  The  stretching  force  constants  are 
higher  than  experimental  values  by  a  factor  of  approximately  two.  The 
orbital  energy  levels  for  ammonia  and  for  oxygen  difluoride  as  a  function 
of  bond  angles  were  in  reasonable  agreement  with  the  full  LCAO-SCF 
Calculations.  A  variation  of  parameters  to  find  a  new  set  of  parameters 
was  attempted  with  unsuccessful  results. 

The  Extended  Hiickel  Theory  (EHT)  can  be  applied  to  aliphatic 
hydrocarbons  but  its  application  to  non-hydrocarbon  molecules  is  quite 
limited • 
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I.  INTRODUCTION 

There  are  many  levels  of  sophistication  for  the  calculation  of 
electronic  structure  using  molecular  orbital  methods.  The  simplest 
approach  is  the  Huckel  Molecular  Orbital  Theory  (HMO).  This  is  one  of 
the  semiempirical  methods  which  is  used  extensively  in  the  quantum 
mechanical  study  of  the  mobile  electrons  of  conjugated  molecules.  The 
HMO  was  extended  by  many  people.  One  of  the  theories  examined  in  this 
paper  is  the  "Extended  Huckel  Theory  (EHT)"  which  was  developed  by 
Roald  Hoffmann  (l).  A  higher  level  of  sophistication  is  the  self- 
consistent  field  (SCF)  molecular  orbital  calculation.  One  of  those 
methods  which  has  recently  been  developed  is  based  on  a  theory  intro- 
duced by  J.  A.  Pople:  "Approximate  Self-Coiisistent  Field  Molecular 
Orbital  Theory  (ASCFT)"  (2),  (3),  (U).  This  ASCFT  Calculation  is 
predominantly  used  throughout  this  text. 

We  can  say  molecular  orbital  methods  are  "curve  fitting".  They 
are  really  based  on  semiempirical  methods.  The  purposes  of  the  calculation 
described  herein  are  the  prediction  of  potential  surface  by  Pople' s 
theory  to  examine  and  to  try  to  improve  the  agreement  of  a  calculated 
potential  surface  with  an  experimental  potential  surface. 

The  background  describing  each  of  these  theories  will  be  discussed 
in  the  next  section. 

II.  Httckel  Molecular  Orbital  Theory  (HMO) 

In  this  theory,  the  mobile  electrons  can  be  treated  as  occupying 
a  set  of  delocalized  molecular  orbitals  (by  the  Pauli  Exclusion  Principle, 
not  more  than  two  electrons  in  each  molecular  orbital).  These  molecular 
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orbitals  are  eigenf unctions  of  a  one  electron  Hamiltonian.  The 
orbitals  are  approximated  by  a  linear  combination  of  atomic  orbitals 
(LCAO)  centered  on  various  atoms.  The  Slater  type  atomic  orbitals 
(£)  are  generally  used  as  approximations  to  the  actual  atomic  orbitals 
for  each  atom.  When  this  theory  is  applied  to  a  wide  range  of  molecules, 
it  has  the  advantage  of  great  simplicity  but  serious  discrepancies  result 
because  of  the  neglect  of  interelectronic  repulsions  for  the  case  of 
many  electrons. 

For  the  unsaturated  molecules :  The  X  -  molecular  orbitals  are 
formed  by  a  LCAO  of  the  2lZ  atomic  orbitals  which,  it  is  assumed,  share 
the  same  nodal  plane.  In  W~X  electron  systems,  the  molecular  orbitals 
are  given  by 

where  LUr  represents  the  j^n  molecular  orbital,  ^represents  the  atomic 
orbital  in  the  ~2J    atom  and  Ct^  represents  the  coefficients  of  the  j^n 
molecular  orbital  and  the  2J       atom  (atomic  orbital).  It  is  noted  that 
the  subscript  i,  j,  k,  1  etc.  stand  for  the  molecular  orbitals  and  ft  j1> ', 
CT  A  etc.  stand  for  the  atoms  or  atomic  orbitals. 

From  a  given  electron  Hamiltonian,  Equation  (l.l)  must  satisfy 
the  equation 

ft%  --  im  £« 

Where  Ej  is  the  eigenvalue  of  the  j"1  state,  \T/-  is  a  column  matrix, 
and^  is  the  Hamiltonian  operator.  This  can  be  put  in  matrix  form: 
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W 


(1.2) 


We  multiply  equation  (1.2)  by  ijj»  and  integrate  over  all  space,  mT  9 
then  we  obtain 


IJj'^dT.jVEjf^T 


(1.3) 


Substituting  equation  (1.1)  into  (1.3) 


Defining 
and 


gCiCj^<MrsgCKA/^ 


H^=/^>4,dr 


We  may  write  equation  (l.U) 


(1.1+) 

(1.5) 
(1.6) 


Taking  the  derivative  of  equation  (1.7)  with  respect  toUfc^  to 
the  energy^   ydLfcX^ J  ,  we  obtain 


(1.7) 


minimize 


It  may  be  shown  that 


(1.8) 

(1.9) 
(1.10) 


Therefore  we  obtain  the  equation 


or 


ICf^(^"-Efe^)  -o 


(1.11) 


(1.12) 
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This  equation  (1.12)  is  the  well-known  secular  equation  and  gives  rise 
to  the  secular  determinants 


I      Hm>  ~~  \&^>m*  I   "   0 


\/KX         Mk^Pil  -  W  (1#13) 

Equation  (1.13)  is  the  non-trivial  solution  of   (1.12). 

Now  we  discuss  the  Hermitian  Property  of  the  Matrices    (  HyuA  ) 
and  \k)jx?\)  •     A  matrix  H  is  said  to  be  Hermitian  if  it  is  equal  to  the 
conjugate  of  its  transpose,   i.e., 

W    =_H  _ 

Where        H=(lV)     <f»i     H'=(Hj>u)  (1-lll) 

By  definition,  a  real  Hermitian  matrix  is  symmetric,  but  an  imaginary 
symmetric  matrix  is  not  necessarily  an  Hermitian  matrix.  It  is  noted 
that  the  eigenvalues  of  an  Hermitian  matrix  are  all  real. 

Next  we  discuss  some  transformation  properties.  In  general,  the 
atomic  orbital  (J)^  does  not  form  a  complete  orthonormal  set.  This 
basis,  (p^  ,  can  be  transformed  into  an  orthonormal  basis  by  using  a 
unitary  transformation  matrix  such  that 


* 


ne*  ■  oU 


i  new 


--ZXA 


*      ~  <^  oftr*  (i.i5) 


where  ty^,    ^s  ^^e  new  orthonormal  basis  formed  by  transforming  the  old 

basis  Cp^  in  N  dimensional  space  and   |2^  is  the  transformation  matrix 

element.  The  old  H  can  also  be  transformed  into  a  new  basis  such  that 

-new        9        +         ( 


hj«»  ;  Z  "pT#J  K#h  dT  o.") 


It  may  be  recalled  that  the  old  basis  is    H^  =  J  ^%  f*t  ty^     ~^ 
Therefore  equation   (1.16)  can  be  rewritten 


\^--2t2W™Tv*  o.i6). 
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h>2/  ~&   '<"*  n«  '^  (1.17) 

It  is  noted  that  if  the  transformations  were  made: 

CJT  -  Z  <feT^  daw 

then  the  equation  (1.17)  becomes 

i*mi>  z  Z_»    Ltr  no  Uz/  (1.19) 

Equations  (1.17)  and  (1.19)  can  be  written  by  matrix  notation  such  that 
new  — | — *r  .  .old.^/ 

H    =T  H  T  (1-17)I 

and 

,  ,  new  -j—*'  i  |  old-,— 

H     :T    n    T  (Li?)- 

From  equation  (1.19)',  if  T  is  a  unitary  matrixf    zr~|   J  then  we  have 

Hlcw=T"lH0ldT  (1-20) 

This  equation  (1.20)  is  frequently  used  and  usually  T  is  real  and  Hold 
is  a  symmetric  matrix. 

Referring  to  the  secular  equation  (1.12),  we  define /Vft#=b-fe*  then 

This  can  be  expressed  in  matrix  notation  in  the  following  manner 

HC-SQ-0  (1-21)1 

c = Co,*) 

Therefore,  we  obtain 

*=OwO  =  C(s"H)C  (1-22) 
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(1.23) 


It  is  generally  true  that  the  matrix  (S  H)  may  not  be  symmetric.  In 
equation  (1.22),  the  eigenvalue  cannot  be  obtained  by  methods  usually 
employed  to  diagonalize  Hermitian  matrices.  One  way  to  solve  this 
problem  is  to  transform  the  matrix  S  to  an  orthonormal  basis,  then 
diagonalize  the  H  matrix.  The  procedures  are  described  as  follows: 
Using  the  transformation  matrix  Q  such  thatllJrQd),  we  transform  the 
matrix  S  to  the  identity  matrix  I 

Sn^I=QS°V 

Also  we  diagonalize  the  matrix  S  by  a  unitary  transformation  matrix  T 
such  that 

(1.25) 

-  [CT^T^T^T)]"'  (1'26) 

Equating  equation  (1.2U)  and  (1.26),  we  obtain 

Then  we  transform  the  old  basis  H  to  the  new  basis  by  the  transformation 
matrix  Q  such  that 

=  (t^t)(h0|JXtW  (1'28) 

Finally  we  can  diagonalize  Hnew  by  the  matrix  U.  Then  we  obtain  the 
eigenvalues 

e=UHneV 


or 


Hi 


where  U  is  the  matrix  which  diagonalizes  the  Hnew  and  T  is  the  unitary- 
matrix  which  diagonalizes  S.  Also  the  correct  eigenf unctions  are 

W  *UTX*T$  a.3o) 

The  general  operations  for  the  application  of  HMO  theory  to 
specific  problems  are 

(a)  Setting  up  the  Secular  determinant 

(b)  Expanding  the  Secular  determinant 

(c)  Finding  the  eigenvalues 

(d)  Determining  the  M.O. 's  by  finding  the  eigenvectors. 

III.  Extended  Hu'ckel  Theory  (EHT) 

Huckel  calculations  originally  were  applied  exclusively  for  the 
Planar  conjugated  and  aromatic  systems  of  J\Z    electrons  with  the  neglect 
of  overlap  for  non-adjacent  atoms.  In  the  extended  Hiickel  theory  (EHT), 
Roald  Hoffmann  (1)  expanded  the  basis  set  to  consist  of  2s  and  2P 
orbitals  for  Carbon  atoms  and  Is  orbitals  for  Hydrogen,  with  the 
inclusion  of  all  overlap  integrals.  The  EHT  can  calculate  molecular 
orbitals  for  the  aliphatic  and  aromatic  compounds  which  gives  good 
results  for  most  hydrocarbons  because  of  the  resonable  set  of  parameters 
chosen  for  the  calculations.  The  parameters  are  nicely  adjusted  to  the 
self-consistent  value  for  the  hydrocarbons  but  when  the  EHT  is  applied 
to  non-hydrocarbons,  it  is  not  as  successful,  and  even  predicts  a  lack 
of  minima  in  many  bond  stretching  motions. 
III-A.  Method  of  Calculation 

For  the  calculation  of  the  molecule  CnHm,  we  use  a  basis  set 
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consisting  of  m-Hydrogen  Slater  orbitals,  n-2s  and  3n-2P  orbitals. 
The  Slater  exponents  were  selected  to  be  1.0  for  Hydrogen  and  1.625 
for  Carbon  throughout  the  calculations.  It  is  noted  that  1.20  was 
selected  for  Hydrogen  in  the  calculation  of  the  ethane  molecule.  Those 
values  are  listed  in  Table  1.  The  order  of  the  matrix  is  Un  +  m  which 
gives  the  number  of  energy  levels  of  the  system  and  the  eigenvalue  E's 
for  off  diagonal  matrix  elements  are  retained  in  the  Secular  Equation. 


TABLE  1 
Slater  exponents  for  valence  Shell  S  and  P  atomic  orbitals. 


H  1.0* 

Li  0.65 

Be  0.975 

B"  1.125 


3 


1.30 


c-  1.U5 

C    1.625 
0+   1.80 
1.95 
2.125 


N 


0" 

2.10 

0 

2.275 

0+ 

2.U5 

F 

2.60 

N 


*  Present  calculations  we  carried  out  using  1.2  for  Hydrogen. 


The  diagonal  elements  <]pM     are  chosen  to  be  valence  state  ionization 
potentials;  the  particular  values  used  are  essentially  those  of  Skinner 
and  Pritchard  (6).  For  example,  the  SP^  valence  state  ionization 
potentials  for  Carbon  and  Is  Hydrogen  atoms  have  the  following  values: 

H^tep) --M.4     e* 

\-\*m(Cx)z-214     ex  (2#1) 

H^  (Ms}  =-13.6     e.y. 

The  choice  of  the  off-diagonal  term    V\mx>    is  approximated     as 
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>ju*  -  0.5  K  (  tyyi  +  Hw)  (S^) 


(2.2) 

The  approximation  (2.2)  was  first  discussed  by  Mulliken  (7)  and  used 
in  the  M.O.  Calculations  by  Wolfsberg  and  Helmholtz  (8). 

The  value  of  K  in  equation  (2.2)  was  chosen  to  be  1.75.  These 
calculations  give  gairly  good  results  for  hydrocarbons. 

IV.  Approximate  Self -Cons is tent  Field  Molecular  Orbital  Theory  (ASCFT) 
Approximate  methods  for  obtaining  self-consistent  field  molecular 
orbitals  have  recently  been  developed  by  J.  A.  Pople  (2),  (3),  (li). 
Slater  type  orbitals  (S>)  are  used  as  basis  functions.  Overlap  integrals 
between  two  different  atoms  were  calculated  using  the  formula  given  by 
Mulliken,  Rieke,  Orloff  and  Orloff  (9).  It  is  assumed  that  the  Slater 
type  atomic  orbitals  are  real  atomic  orbitals,  and  that  the  M.O. 's  for 
the  ground  state  are  obtained  by  Linear  Combinations  of  Slater  type 
atomic  orbitals.  If  the  Coefficients  in  the  LCAO  Molecular  orbitals 
are  chosen  to  minimize  the  total  energy,  one  obtains  LCAO  -  SCF  orbitals. 
The  best  LCAO  approximation  to  SCF  functions  were  first  given  by  Roothaan 
(10).  It  is  noted  that  the  EHT  treated  all  the  electrons  independently. 
EHT  wave  functions  are  used  as  a  starting  point  for  the  self -consistent 
field  interations  in  this  ASCFT. 

The  SCF  methods  were  developed  using  two  approximations  based  on 
the  neglect  of  the  differential  overlap  between  all  valence  atomic 
orbitals.  One  of  these  involves  complete  neglect  of  differential  overlap 
(CNDO)  in  all  basis  sets  connected  by  the  transformations  as  mentioned 
in  Section  I,  that  is :  the  of f -diagonal  elements  of  overlap  matrix 
in  the  Hamiltonian  matrix  are  to  be  neglected.  The  other  involves  the 
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neglect  of  the  products  of  (4L0)  0^.0)  °^  any  ^w0  atomic  orbitals 
(D^   and  (h       in  all  electron-repulsion  integrals.  This  complete 
neglect  of  differential  overlap  of  all  basis  sets  will  result  in  non- 
orthonormal  basis  sets. 

The  CNDO  method  in  ASCFT  should  be  regarded  as  an  approximation 
to  a  full  LCAO  -  SCF  calculation  which  uses  a  minimal  basis  containing 
only  inner-shell  and  valence-shell  atomic  orbitals.  The  Slater 
exponents  used  in  these  calculations  are  listed  in  Table  1  (for  Hydrogen, 
1.2  is  used).  The  values  f or  /O   and  ionization  potentials  in  2s  and 
2p  atomic  orbitals  from  Lithium  to  Fluorine  were  set  semiempirically. 
Those  values  used  in  these  calculations  are  tabulated  in  Table  2  and 
Table  3  respectively. 


TABLE  2 

Bonding  Parameters    /J a     (in  e.v.) 

Atom       H       Li       Be        B        C          NO 

F 

9       9        13      17       21      25       31 

39 

TABLE  3 

"Average"  Ionization  Potentials  A.n  (in.  e.v.) 

Atom     H       Li     Be      B       C  NO 
Is     13.06 

2s     —     5.39    9.32    1U.05    19.  Ui  25.58   32.38 

2P     —     3.5U    5.96          8.30    10.67  13.19   15.85 


U0.20 
18.66 
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The  CNDO  method  of  ASCFT  can  be  used  for  closed  ground  shell 
structures  and  doublets,  but  a  triplet  structure  is  not  available  with 
the  computer  program  used  (for  example  O2,  Bg  and  C2)»  The  computer 
program  is  presented  in  Appendix  1. 
IV-1.  Self -Consistent  Field  (SCF)  Equation 

The  initial  discussion  covers  Roothann's  LCAO  -  SCF  method  (10) 
as  a  basis  for  calculating  the  MO's  of  a  molecule,  paying  particular 
attention  to  transformation  properties.  The  Is  Hydrogen  and  valence 
shell  electrons  for  Li  to  F  are  assigned  to  LCAO  molecular  orbitals. 
Minimization  of  total  energy  by  a  variation  principle  with  SCF  orbitals 
leads  to  a  series  of  simultaneous  non-Linear  equations,  the  so  called 
Hartree  Fock  equations  (11),  which  are  applied  to  LCAO  molecular 
orbitals  with  a  closed  shell  configuration.  These  equations  were  shown 
by  Roothaan  to  reduce  to  the  form  (10) 


2-  v/j&LvX    -  J?  ^#G//t£/c 


(U.D 


where  (using  atonic  units) 

VfiD  -  ykv  +G/nv  (u.2) 


0*4$  - J  j  j)P  4^  "fir  k<*  fa21  **  fe        <*.* ) 


(U.3) 
(U.U) 
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S>ua    =  jfy?  9*  to 


(U.6) 


occ 


A. 

Where  Ei  is  the  orbital  energy  for  the  i™1  M.O.  VyXV   is  the  matrix 
element  of  the  one  electron  Hamiltonian  including  the  kinetic  energy 
and  the  Potential  energy  in  the  electrostatic  field  of  the  core,  this 
being  written  as  a  sum  of  potentials  V^(r)  for  various  atoms  A  in  the 
molecule.  vTui»  is  the  matrix  element  of  the  potential  due  to  valence 
electrons  and  depends  on  the  M.O. 's  via  the  population  matrix  nac"    . 

The  total  energy  of  the  molecule  is  obtained  by  adding  the  repulsion 
energy  between  isolated  cores  and  the  separated  valence  electrons.  The 
repulsion  term  can  be  approximated  by  a  point  charge  model 

Etotal    =  E electron;^ EZ  *$*M*  (U.8) 

/KB  / 

where     2_jR      is  the  core  charge   of  atom  A.     The  electronic  energy  is 


given  by 


A'J  (U.9) 


/*.*         '  u  xt 
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(U.io) 

17-2.  Theory  with  Complete  Neglect  of  Differential  Overlap  (CNDO) 

First  consider  the  simplest  version  of  the  theory  which  involves 
the  complete  neglect  of  differential  overlap  between  atomic  orbitals 
on  the  same  atom.  The  approximations  are  given  below. 

Approximation  A 
The  atomic  orbitals  (Jpn     are  treated  as  if  they  form  an  ortho- 
normal  set  such  that  Oi|2/  *  CU2f      an<^  then  the   L.>uU    forms  an 
orthonormal  matrix.  The  condition  for  orthonormality  becomes: 


lC^m  =  $ij 


j*        -     *  u  Oull) 

With  this  approximation,  the  diagonal  matrix  elements  \]mju         (electron 
density)  correspond  to  the  electron  populations  of  the  atomic  orbitals 
%      and 

where  2N  is  the  number  of  valence  electrons. 

Approximation  B. 
Equations  (U.U)  are  neglected  because  of  the  orthonormality  of 
atomic  orbitals.  This  means  { /\^J\7^   )  is  zero  unless  JL\zz2J      and 
/\=S~  .  The  non-zero  values  will  be  written   (M/if  such  that 
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^6?yu    =  (**(>"/■) 


(U.13) 


Approximation  C 
The  electron- interaction  integrals  (UiZ;  are  assumed  to  depend 
only  on  the  atoms  to  which  the  orbitals  (jV,  and  (J^,  belong  and  not 
on  the  actual  type  of  the  orbitals.  This  means  that  there  remains  only 
a  set  of  atomic  electron-interaction  integrals  OAb     measuring  an 
average  repulsion  between  an  electron  in  valence  atomic  orbital  on  atom 
A  and  another  in  a  valence  orbital  on  atom  B.  The  set  of  atomic  orbitals 
on  one  particular  atom  A  is  replaced  by  an  alternative  set  tm3 

X     r  ™  Ouiu) 

A  general  electron-replusion  integral  in  this  new  basis  set  is  given  by 
Now  we  apply  approximation  C,  then  equation  (U.l£)  becomes 


(  tm+» [  <$>*  (j)J  =  Ya&  Smr\  (U.16) 


It  is  noted  that  equation  (I4..I6)  is  invariant  under  the  orthogonal 
transformation  matrix  0  which  gives  the  new  basis  x». 

Using  approximation  A,  and  B,  the  matrix  elements  Y/AlJ    become 


*>a 


-      hj^  "~k~\/iM  \jm  +Zi  f^cr^fi" 


(U.17) 
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hi*      =  fy*  "Tljw"^* 


Equation  (U»17)  and  (U»l8)  can  be  rewritten  by  approximation  C 


l>  =    ^2 6^^a+Ba4+Z  Mb  (U9) 


f>^    z     n/a,  T  p    Rm^B 


(U.20) 


where  yM  belongs  to  atom  A  and  Pgg  is  the  total  valence  electron 
density  on  atom  B. 


IBB  ~  /  f    ri. 

Also  the  core  matrix   H>u>|  can  b®  separated  as  follows: 


(U.21) 


B#0 


=  LU  -£  C^IVbW 


there      U/uj       is  the  diagonal  element  of    Cpw     with  respect  to  the 
one-electron  Hamiltonian  containing  only 

the  core  of  its  own  atom  (  LJ^u/u  is  an  atomic  quantity  measuring  the 
energy  of  the  atomic  orbital)  and  VJLoj  is  the  one-electron  matrix 
element  using  the  local  core  Hamiltonian  (  JA     and   ~2J      are  on  the 
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same  atoms).  It  is  noted  that  the  cores  consist  of  the  nucleus  and 
inner-shell  electron. 

Approximation  D 
The  interaction  terms  of  the  distribution  hju  Q^>   on  atom  A 
and  the  nucleus  of  atom  B,  (  jU  Vf3  ^  )  are  equal  to  zero  if  /J  ZE-  2J  . 
The  integrals  (M  |  Vfc  |  ^  )  are  taken  to  be  the  same  for  all  valence 
atomic  orbitals  on  Atom  A  and  there  results 

VjiB   =  \P    \k\Aj  (U.2U) 

It  should  be  noted  that  the  attraction  term  matrix  element  v^b 
between  nucleus  B  and  the  adjacent  valence  electron  M       in  the  field 
of  the  nucleus  A  need  not  be  symmetric.  Therefore  the  equations  (1+.22) 
and  (U.23)  can  be  rewritten 

\~\MM      =  UmM    T  2-J   V/4B  (  U  on   *4e»  A   J 

It  is  noted  that  rj/i^^-0  if    <Pm       an<^  ^^   are  on  different 
atoms  which  will  be  discussed  in  the  next  approximation. 

Approximation  E 
If   CD^  and  (p2)  are  on  different  atoms,  it  is  only  necessary  to 
consider  the  possible  lowering  of  energy  levels  by  being  in  the  electro- 
static field  of  two  atoms  simultaneously  which  is  referred  to  as  a 
"Resonance  Integral"  and  is  denoted  by  /jm~v     .  \yA~D       can  be 
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approximated  as 

f>2  =r  fyn>    */£&  ^jaV  (U#27) 

W'lere  /0/*b   is  a  parameter  depending  only  on  the  nature  of  the  atoms 
A  and  B.  Empirical  values  will  be  used  for  these  parameters. 

The  Yuzu         matrix  elements  now  reduce  to  the  form  from  equation 
(U.19)  to  (U.27) 

^  B(*A)  ^  (U.28) 

where  (|X,  belongs  to  atom  A  and  (f^,  belongs  to  atom  B. 

Using  the  same  approximations,  the  total  energy  is  given  by  one 
or  two  atom  types  from  the  equation  (I4..8) 

L,  MA  =LEa   +  ZZ  &6 

Where 


B 


(U.30) 
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IV-3.     Choice  of  Parameters  of  CMD01  and  CND02 

In  the  complete  neglect  of  differential  overlap  approximation, 
the  CND01  and  CND02  methods  are  different  levels  of  approximation 
of  the  equations. 

The  LCAOSCF  Hatree-Fock  matrix  elements  are 


Bfc4) 


(a.28) 


where  the  bonding  parameter    yOAQ      is  given  by 

j3°B  =  6.5  (/3l+ySl^>  a.32) 

In  the  equation  (U.32),  /Sa    ^s  chosen  empirically  but  depends  only 
on  the  nature  of  atom  A.  These  values  are  listed  in  Table  2,  The 
parameters  L-/2S,2  5  and  U^P32p  can  be  estimated  from  ionization 
potentials,  using  the  theoretical  value  of   OXX  (12). 

or 

and  the  atomic  electron  affinities        /\n  S    are  used  by  the  formula 

f\/X    =     U/4/A      +  JjAOAA  (U.36) 
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Also  the  equation  (U.28)  can  be  rewritten 


V/au   ~  Dam  +(P/irTW  ~%AA 


y/4/4 


BOM) 

where  the  last  term  of  equation  (U.28)'  is  described  as  penetration 
integral  contributions  to  rjU/U       . 

The  following  table  gives  the  approximate  form  between  CND01  and 
CND02. 


reference 

-U 


TABLE  k 
Difference  between  CND01  and  CND02 
CND01  CND02 


'/*/* 


Ou37) 

P(*)l   ^°      '       B^>         (U.38) 
Where  1  the  average  ionization  potentials  Xm     (©•▼•)  used  to  fix 


r>w 


>tyM 


u. 


'/»/» 


are  listed  in  Table  3. 


2  the  matrix  elements  J(ls  +  As)  and  ^(ip  +  Ap)  are  listed 
in  Table  5. 


TABLE  5 

Matrix  Elements  From  Atomic  Data 

(e.v.) 

H             Li               Be              B                C 

N 

0              F 

i(la+Aa) 

7.176      3.106      5.9U6      9.59U      lU.Ctfl 

19.316 

25.390     32.272 

i<W 

1.258       2.563      U.001        5.572 

7.275 

9.111    11.080 
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If  the  penetration  integral  term  in  equation  (U.28)»  was  neglected 
to  avoid  the  excess  bonding  and  high  value  of  force  constants,  then 
equation  (U.28)'  can  be  reduced  to  the  form  of  (U.38).  This  corresponds 
to  substituting  the  following  term  for  CND02  from  CND01: 


Vab  =  2bTa&  (^39> 


Using  parameters  specified  in  the  above  way,  the  LCAOSCF  equations  are 
solved  in  the  following  steps; 

(a)  An  initial  set  of  M.O.  coefficients  (_/LyU   is  obtained  from 
a"Huckel  type"  calculation  such  that  r/j/U   is  replaced  by 
average  ionization  potential  and  Pajoj       by  A  SmU  * 

(b)  Electrons  are  assigned  in  pairs  to  the  M.O. 's  with  lowest  energies. 

(c)  Calculated  \/yu     and  Y/jci)       are  used  to  form  a  new  Hatree  Fock 
matrix  element  YJA2J         • 

(d)  New  coefficients  (.Lys  are  derived  from  this  rjjn>       and  the 
process  is  repeated  from  step  (b)  until  self-consistency  is 
achieved  on  all  coefficients  with  a  tolerance  of  0.0001. 

V.  Discussion  of  Results 

The  molecules  selected  for  these  calculations  are  as  follows: 
Hydrogen  molecule  (Hr>),  Lithium  molecule  (Li2),  Beryllium  (Bez), 
Nitrogen  molecule  (Nz),  Fluorine  molecule  (Fz),  water  molecule  (J^O), 
Ammonia  molecule  (NH3),  Ethane  molecule  (C2E^) }     oxygen  difluoride 
(F20),  Difluorodiazine  (N2F2),  Methane  (CH^),  Methylfluoride  (CH3F) 
and  Hydrogen  Peroxide  (l^O^)*  These  calculations  will  show  the  minimum 
configurations  for  bond  lengths,  bond  angles,  force  constants,  the 
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energy  level  considerations  for  NH3  and  F2O,  and  the  Barrier  for  NH3, 
C2H6,  N2F2  and  H2O2.  Finally  the  calculations  will  extend  the  variation 
of  parameters  to  fit  the  experimental  potential  surfaces. 
V-l.  Comparison  of  the  equilibrium  configurations  and  force  constants 

The  ASCFT  calculations  with  various  inter-nuclear  distances  and 
bond  angles  are  useful  in  predicting  equilibrium  bond  lengths  and  bond 
angles.  These  results  are  listed  in  Table  6.  The  equilibrium  configur- 
ations (bond  lengths  and  bond  angles)  are  in  fairly  good  agreement  with 
experimental  values  but  the  force  constants  are  too  high  compared  with 
the  experimental  v alues  by  a  factor  of  about  2.0.  This  can  be  improved 
by  adjusting  parameters,  such  as  the  Slater  exponent  and  yO        . 
Sample  results  are  discussed  in  Section  V-9. 


Comparison  of 

TABLE  6 
Bond  Lengths,   Bond  Angles  & 

Force  Constants   ("Rf  «  "ft®) 

(13) 

Bond  Length  (A0) 

Bond  Angle   (cleg) 

(13^ Force 
Consts 

Ue:  Id"  erg /ma?  J 

Cal. 

Obs. 

Cal. 

Obs. 

Cal. 

Obs. 

N2 

0.7U59 

0.7U15 

10.11*3 

^.6 

1.1395 

1.0976 

52.160 

23.10 

H20 

1.0289 

0.958 

10U.22 

10lu  k$ 

17.260 

8.23     ,-,  v 

0.81*5 

0.69U  {1U) 

NH3 

1.0657 

1.015 

lll.68 

6.86     .  , . 

106.825 

106.6 

0.95 

0.60     C11*) 

0% 

1.1135 

1.1135 

12.65 

**U5     fin 

109.U7 

109.1*7 

O.U98  (15) 

Rcc 
J  iccw 

1.1195 

1.102 

12.12 

5.35 

1.U633 

1.5ii3 

16.31 

U.57 

107.76 

109737' 

0.51    l*> 

110.86 

109°37' 

___________ 

0.63 

-      — < 
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V-2.  Water  molecules  (H2O) 

For  H2O,  C2U  symmetry  was  ass  Timed  with  an  0-H  bond  length  of 
O.96A0.  The  first  calculation  shows  an  energy  minimum  at  a  bond  angle 
of  107.28°.  At  this  bond  angle ,  the  equilibrium  0-H  bond  length  of 
1.0289A0  is  obtained.  With  this  bond  lengthy  the  second  calculation 
shows  the  bond  angle  of  10h.22°  at  which  the  0-H  bond  length  is 
calculated  to  be  1.0289A0.  This  is  in  satisfactory  agreement  with  the 
experimental  bond  angle  of  10U.li5o  (13 )•  These  calculations  demonstrate 
the  0-H  bond  stretching  and  angle  HOH  bending  motion  such  that  the  0-H 
bond  lengths  tend  to  have  symmetrical  stretching  motion  and  the  bond 
angle  tends  to  have  scissoring  motion  (symmetrical  in-plane  deformation), 
A  series  of  these  calculations  is  listed  in  Table  7.  The  stretching 
force  constant  is  calculated  to  be  17.26  x  10-5  (  ^rne/cm)  which  is 
higher  by  a  factor  of  2.1  with  respect  to  the  experimental  value  of 
8.23  x  lO^  (dyne /cm).  The  bending  force  constant  is  calculated  to  be 
0.8U5  x  10-11  (erg/rad2)  which  is  higher  by  a  factor  of  1.22  than  the 
experimental  value  of  0.69k  x  lO"11  (erS/rad2)  (i^). 


TABLE  7 

Calculation  of 

RoH 

and             HOH  for  H20 

Calculation  I 

Calculation  II 

Experiment 

*oH(A°) 

0.96   (Exp) 

1.0280 

0.958 

ZHOH 

107.28 

10U.22 

10U.U5 

*  The  computer  outputs   (sample) 

are  listed  in  Appendix 

2. 

V-3.  Ammonia  Molecule  (NH3) 

The  ammonia  molecule  with  an  NH  bond  length  of  1.05A°  (13)  is 
calculated  to  be  non-planar  with  an   ^  HNH  bond  angle  of  106.825° 
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which  is  in  close  agreement  with  the  experimental  value  of  106.6°  (13)« 
At  this  value  of  the  bond  angle,  the  NH  bond  length  was  found  to  be 
1.0657  A°.  The  subsequent  calculation  shows  a  bond  angle  of  10U.860 
and  again  the  bond  length  of  1.0673  A0.  On  the  other  hand  the  equili- 
brium bond  length  for  the  planar  configuration  gives  1.05U3  A°.  It  is 
interesting  when  calculations  presented  here  are  compared  with  results 
recently  calculated  by  R.  Moccia  and  L.  Randaccio  (17).  They  reported 
the  minimum  bond  length  of  I.878  a.u.  (0.9938A0)  at  a  bond  angle  of 
105°36«  and  a  bond  length  of  1.837  a.u.  (0.9721  A°)  at  a  planar 
configuration.  The  change  in  minimum  bond  length  between  non-planar 
and  planar  configuration  is  0.0217  A0  by  Moccia  and  0.011U  A°  by  the 
ASCFT.  The  smaller  value  of  this  bond  length  change  in  the  present 
calculation  is  due  to  the  higher  stretching  force  constants  for  ASCFT. 
Also,  subsequent  calculations  can  predict  the  symmetrical  bond  stretching 
motion  and  symmetrical  in-plane  deformation  of  bond  angle  for  the  ground 
state  of  this  molecule.  The  stretching  force  constants  are  calculated 
to  be  lh.6h  x  10^  (dyne/cm).  This  differs  by  a  factor  of  2.1U  with 
respect  to  the  experimental  value  of  6.86  x  10^  d7ne/cm  (Hi)  and  the 
bending  force  constants  are  calculated  to  be  0.95  x  10"*11  erg/rad.2 
(factor  of  1.58  with  respect  to  experimental  values  of  0.6  x  10" 
erg/rad2  (Hi)). 

It  is  also  interesting  to  compare  the  barriers  calculated  by  the 
several  methods  and  the  difference  in  barriers  with  respect  to  adjusting 
the  various  parameters.  These  are  listed  in  Tables  8  and  9,  respectively. 
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TABLE  8 
Comparison  of  the  results  of  Calculations  for  NH3 

Reference 

Method  of  calculations 

RNH  fo.ll) 

^HNH(dq) 

Pyramid           Barrier 
Height(a„u.  )Height(a.u.> 

Present  Calc, 

Approximate  (CNDO) 
SCF   MO 

1.9923 

106 °U9» 

0.7U9 

0.0251 

Higuchi  (18) 

Approximate  polycenter 
SCF   MO 

1.916 

108.0 

0.592 

0.01U7 

Kaplan  (19) 

Single  Determinant 
polycenter  SCF  MO 

1.916 

106°^, 

0.719 

0.1U6 

Moccia  (20) 

Single  Determinant 
one  center  SCF  MO 

1.928 

110°2 « 

0.62H5 

0.0027 

Joshi   (21) 

Single  determinant 
one  center  SCF  MO 

1.867 

109°3U» 

0.626 

O.OUll 

Moccia  (17) 

1.878 

10£°36» 

0.739 

0.00905 

Experiment 

i 

1.9162 

106 °U7' 

0.719 

0.0092 (22) 

From  Table  9,  it  is  interesting  to  compare  the  change  in  barrier 
with  respect  to  a  change  in  bond  angles  at  the  fixed  bond  length  of 
1.0657A  (Calculation  II  and  IV)  such  that  the  change  in  barrier  is 
0.000030279  a.u.  (0.0025U0830  -  0.00251055D  and  the  change  in  bond 
angle  is  1.965°  (106.825°  -  10U.860).  From  this  calculation,  it  is 
found  that  the  calculated  barrier  is  0.0025U0830  a.u.  which  is  different 
by  only  1,2$  from  the  first  calculated  barrier. 

The  electronic  configuration  is  found  to  be 

(lai)2  (le)**  (2a1)2  ^  # 
Configurations  calculated  by  ASCFT,  EHT  and  full  LCAOSCFT  by  H.  D. 
Joshi  (21)  are  compared  in  Fig.  1  and  Fig  2.  There  is  complete  agreement 
between  ASCFT  and  full  LCAOSCFT  calculations  for  the  orbital  energies 
as  a  function  of  height  of  pyramid.  It  is  noted  that  the  lowest  orbital 
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fc.2.  Ek&G»oF  NH3 vs.  PiWttllD HtlfrHT 
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0.2  0.3  0.4 

HEI&HT  OF  PYRAMID  (A*) 


3h 


energy  in  LCAOSCFT  calculation  is  due  to  the  (is)  electron  (inner- 
shell)  of  Nitrogen  but  ASCFT  neglects  that  part.  On  the  other  hand, 
the  EHT  calculation  (Fig.  2)  does  not  even  demonstrate  a  double  minimum 
potential  curve.  Therefore  the  EHT  is  quite  inferior  for  this  kind  of 
molecule. 


Change 

in  Barrier  for 

RnhC/H 

TABLE 
NH3  w.r. 

9 

t.  minimizing  the  parameter 

2HHW          E total  (tf-U.)        Barrier  (a.u,) 

Calculation  I 
(General) 

1.02   (fixed) 
1.02   (fixed) 

0 
112.0 

120 
106.825 

-13.85761125 
-13.87775028 

0.02013903 

■ 

Calculation  II 
(Equilibrium) 

1.05U3 
1.0657 

0 
112.0 

120 
106.825 

-13.86U30636 
-13.889U12U5 

0.02510551 

III 

1.05U3 
1.05U3 

0 

113.18 

120 
105.10 

-13.86U30636 
-13.8889015U 

0.02U53518 

17 

1.0657 
1.0657 

113?55 

120 
lOlt.86 

-13.86362706 
-I3.8897H166 

0.02608760 

Experimental 

112.206 

106*6 

0*0092 

_.,  The, true  Barrier  for  this  calculation                                             ~  ^^  -0,n 
*  is  (I3.889711i66  -  I3.86U30636)   (a.u.)                                            0.025U0830 

V-U.  Ethane  molecule  (CpHz) 

Calculations  based  on  ASCFT  and  EHT  have  been  carried  out  for 
Ethane  in  both  eclipsed  and  staggered  configurations  with  experimental 
values  such  that  ^  =  1.5U3A0,  Bf£     =  1.102A0  and  J_  CC  H   =  109°37« 

(tetrahedral  angle).  The  first  calculation  shows  the  equilibrium  angle 
for  both  eclipsed  and  staggered  forms  using  EHT  and  ASCFT.  The  energy- 
then  was  minimized  by  varying  the  bond  length  R^£,  and  finally  by  varying 
bond  length  R^j.  The  stretching  force  constants  ( ~KCC   and  KCH  )  were 
found  and  those  values  are  listed  in  Table  10. 
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It  is  of  interest  to  look  at  the  barrier  resulting  in  the  change 
in  geometry.  These  changes  in  parameters  and  the  corresponding  barriers 
are  compared  in  Table  11.  These  comparisons  show  that  the  eclipsed  form 
is  more  open  than  the  staggered  form  by  a  few  tenths  of  a  degreej  also 
the  bond  length  R~  is  elongated  a  few  thousandths  of  an  Angstrom  and 
the  bond  length  RcH  is  shortened  a  few  ten  thousandths  of  an  Angstrom. 
Some  researchers  have  shown  an  elongation  for  R^.  Although  these  changes 

(RcB)  in  geometry  are  small,  their  effect  on  the  internal-rotation 
barrier  is  very  significant. 


TABLE  10 

Comparison  of  Ethane 

Betn.  EHT  & 

ASCFT  Calculation 

Reference 

Calculated 

Expt.   (16) 

ASCFT 

EHT 

Eclipsed 

Stagger 

Eclipsed       Stagger 

D3d 

J-CCH 

111.28 

110.86 

110.99 

110. kk 

109°37» 

^MCH 

107.60 

108. 0U 

107.90 

108. U8 

109°37f 

Rcc 

1.U660 

1.1*633 

1.8070 

1.8000 

1.5U3 

RoT 

1.1193 

1.1195 

0.9086 

0.9088 

1.102 

-feiHCH 

(1.7U) 
(         ) 

(1.67) 
(         ) 

(8.85) 
(         ) 

(8.8) 
(         ) 

0.63 
0.51 

-ftcc  r.rf     1 

15.887 

16.3UO 

197.0 

18U.0 

U.57 

AcH  * 

drM 

12.125 

12.121* 

l6h.05 

186.56 

5.35 

Barrier 

2.3288 

2 

,6586 

2. 875^.125 
(23) 

There  is  some  uncertainty  in  the  Table  11  for  EHT  calculation. 
Every  reported  calculation  scheme  has  the  same  tendency  to  increase 
the  barrier  when  the  parameters  are  changed  in  sequence,  but  EHT 
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calculation  shows  the  barrier  is  decreased.  The  resultant  barrier  by 
EHT  is  1.6358  kcaVmole.  This  is  not  in  agreement  with  the  value 
2.6586  c  /mole  found  by  varying  only  one  parameter,  ZCCW        ,   holding 
the  other  two  parameters  R~  and  R~  constant.   It  is  noted  that 
R.  Hoffmann  (l)  reported  U.O  kcal/mole  with  a  Slater  exponent  (H)  of 
1.0.  Present  calculations  use  the  value  of  1.2.  However  the  results 
are  in  fairly  good  agreement  with  experimental  values  using  both  ASCFT 
and  EHT  (with  some  uncertainty). 


TABLE  11 
Comparison  of  Barrier  for  Ethane 

Reference 

/\  Parameter (Eel. — Stag) 

Barrier (Eclip. — Stag ) 

J-cch  (de<i) 

Kcc   (A') 

KCH  (A°) 

ZCCH 

KCC 

RCH 

T.L.Hill^ 
(191*8) 

0.3U 

0.00514 

0.0002 

0.613 

0.665 

0.700 

L.S.  Bar tell tf#) 
(I960) 

0.20 

0.0029 

0.0001 

0.367 

0.383 

0.39U 

J.T.Vanderslice^' 
&  E.A.Mas  on  (I960) 

0.17 

0.002U 

-0.0001 

0.387 

0.1+16 

O.U3h 

L.Salem^2?) 
(1962) 

0.05 

0.0007 

0.0000 

0.10U 

0.106 

0.107 

i ■ 

1.298 

1 

K.E.Hcnaett^28) 
(I960) 

0.U0 

0.0059 

-0.0002 

0.779 

1.108 

(29) 
V.Magnascov   yj 

(1962) 

0.55 

0.0080 

-0.0011 

-1.532 

1.675 

2.9UU 

O.J.Sovers  &  (30) 
M.Karplus   (1966) 

0.20 

0.0029 

2.311 

3.611 

Present 
(ASCFT) 

(EHT) 

0.U2 

0.55 

1 

0.0027 
0.0070 

-0.0002 
-0.0002 

1.8637 
2.6586 

2.28U8 
1.57U6 

2.3288 
1.6358 

Expt.                           3.03±0.3(kcal/mole)   (3D               2.875!0.125  (kcal/mole) 

(23) 

37 


V-5.  Methane  (CH^)  and  Methylfluoride  (CH3F) 

Calculations  between  methane  and  methylfluoride  were  carried  out 
using  tetrahedral  symmetry  for  methane  and  C3v  symmetry  for  methyl- 
fluoride with  R£h  ■  1.1135A0  for  CH^  and  R££  =  1.109  and  R~  =  1.385 
for  CKoF.  The  calculated  values  for  the  equilibrium  configurations  are 
listed  in  Table  12 .  The  force  constants  are  calculated  for  each  molecule 
and  these  are  listed  in  Table  12.  The  stretching  force  constants  are 
different  by  a  factor  of  2.3  for  CH^  and  2.5  for  CHUF  with  respect  to 
their  observed  values  (15). 

For  a  methane  molecule,  the  calculated  bond  lengths  and  bond  angles 
are  the  same  as  the  observed  values. 

For  methylfluoride  molecules,  the  bond  angle  Z  HCH  is  more  open 
than  the  angle  of  <£  HCF  which  can  be  explained  as  follows  s  there  exists 
rehybridization  for  the  CHoF  molecule  because  there  must  be  more 
P-character  in  C-F  bond  due  to  the  electron  withdrawing  fluorine;  on 
the  other  hand,  there  must  be  more  S-character  in  C-H  bond  which  gives 
more  SP^  hybrid  character  for  the  C-H  sigma  ( <y)  bonds.  It  is  obvious 
that  the  SP^  hybrid  orbital  is  more  planar  than  SP3.  Therefore  ^  HCH 
is  more  open.  It  is  here  that  the  best  agreement  between  theoretical 
and  calculated  results  was  found. 
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TABLE  12 
Comparison  of  f  ,    0   and  kr  Bk.  k©  between  CH^  &  CH^F 

CHU 

Bond  Length  (A0)     Bond  Angle  (deg) 

Force  Consts. 

Cal. 

Obs. 

Cal. 

Obs. 

Cal. 

Obs. 

1.1135 

1.1135 

12.65 

$.h$ 

109. U7 

109.U7 

(.7U5) 

O.U98 

CH3F.NCH 

1.1180 

1.109 

12.09 

ii.7l(l5) 

1.3U58 

1.385 

23.57 

5.6U(15) 

110.382 

109.997 

O.UU3C15) 

108.5U 

108. 9k 

(1.77) 

0.792(15) 

V-6.  Difluorodiazine  (N2F2) 

Calculations  have  been  carried  out  for  difluorodiazine  in  both 
trans  and  Cis  configurations  with  Knn  »  1.25A°,  ^nr  =  l.UUA°  and 
-^NNF   =  115°  (13).  In  both  configurations,  the  new  equilibrium 
distances  and  bond  angles  are  determined  and  are  listed  in  Table  ll*. 
When  comparisons  are  made  of  the  equilibrium  geometry  in  Table  ll*,  the 
bond  angle  ^.NNp   is  more  open  for  Cis  than  trans  configurations. 
Significantly,  the  bond  length  Rjjjj  is  elongated  a  few  thousanths  more 
for  trans  than  Cis  (same  for  K  ^F   ). 

It  is  of  interest  to  compare  the  energy  level  difference  between 
Cis  and  trans  configurations.  The  experimental  results  giving  the 
stability  between  both  configurations  are  listed  in  Table  13.  In 
1963>  Armstrong  and  Marantz  published  that  the  Cis  form  is  more  stable 
by  3.0  2  0.3  kcal/mole  (36)  (refer  also  (37)).  The  present  calculation 
predicts  the  Cis  form  to  be  more  stable  by  3.21*9  ^c  /mole  than  the 
trans  (C2h)  form.  Fortunately,  these  results  are  in  most  satisfactory 
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agreement  with  experimental  value.  In  Table  H4.  are  compared  the  diff- 
erences in  energies  between  the  Cis  and  trans  configurations  for  a 
variation  in  molecular  parameters.  First  the  angle  ZNNjh  was  varied 
to  minimize  the  energy.  With  this  bond  angle,  the  minimum  bond  length 
}^N  was  calculated  and  finally  the  minimum  bond  length  Rjjf  was  found.  In 
each  case,  the  difference  in  energies  was  compared  with  the  one  obtained 
initially.  It  was  found  that  there  is  no  effect  due  to  varying  the 
angle  ^NNJT  and  bond  length  KNN  but  there  is  a  very  significant 
difference  when  the  bond  length  r\Np  was  minimized.  This  result  shows 
the  configuration  is  affected  significantly  due  to  the  KNp  stretching 
motion. 


TABLE  13 
Results  of  experiments  for  N2F2 (stability) 


Year 
19k7 


Author 
S.H.  Bauer (32) 


1959 


I960 


C.B.  Colburn,et.al. 
R.H.  Sanborn  ^t~ 


n^ 


1961 


R.  Ettinger,et.al, 


t#t 


1963 


G.T.  Armstrong 
&  S.  Marantz 


T35T 


Present 


Reference 
J,Am.Ghem.Soc.,69,310U 


U   tt    II 


"  ,81,205 


j.Chem.Phys.,  33,1855 


11  ti 


,  3U,2187 


J.Chem.Phys.,  38,169 


Results 

trans  is  more  stable 
than  Cis  configuration 


Cis  and  N=N<p   is 
more  chemically  active 
is  omer 


Cis  is  more  stable  by 
3.0±0.3  kcal/mole 


Cis  is  more  stable 
by  3.21*9  kcalAiole 


The  next  interesting  calculation  is  the  barrier  due  to  the  rotational 
motion  which  was  reported  by  the  Binenboyn  (37)  .  It  was  reported  by 
experiment  to  be  32.0  kcal/mole  and  calculated  to  be  8U  kcal/mole  (36) 
between  trans  and  linear  configuration.  However,  the  present  calcu- 
lations were  carried  out  in  various  ways  and  those  results  are  listed 
in  Table  15.  The  calculated  barrier  between  linear  and  trans  form  in 


UO 


TABLE  li| 
Equilibrium  Configuration  Btn.  G>ir  &  Czfi  for  N0F2 

Calculation 

Parameters 

A  Par  (&r-Gff  )               Barrier11 

^NNp^) 

Rnh  (A0) 

Rhfm°; 

^NNR  J 

Km 

Rnf 

/^(Cv-Czli) 

^enem  1 

115 

115 

1.25 
1.25 

1.1* 
1.1*1* 

0 

0 

0 

0.2553(0) 

II            Gu- 
Gh 

109.98 
101.09 

1.25 
1.25 

l.Wi 

1.25 

5.89 

0 

0 

o.2553(0) 

III       Caur 

Gk 

109.98 
IOU.09 

1.2197 
1.2218 

1.25 
1.2^ 

5.89 

-0.0021 

0 

(-)0.0^6231) 

IV             G^ 

109.98 
lOh.09 

1.2197 
1.2218 

1.2U10 
1.21*1*1 

5.89 

-0.0021 

-0.0031 

(-)3.2li92 
(3.5) 

Expt. 
value 

— — 

(-)3.0i0.3 

1 1. 

Table  15  is  too  high  with  respect  to  the  reported  barrier  but  the 
calculation  between  trans  and  Cfc  form  shows  a  satisfactory  agreement 
with  the  calculated  (38)  and  experimental  value  (39).  It  must  be 
noted  that  when  an  attempt  was  made  to  calculate  the  rotational  barrier, 
no  dihedral  angle  exists j  that  is,  when  the  dihedral  angle  is  90  degrees, 
the  bond  angleZNNp  turns  out  linear,  and  also  the  dihedral  angle  does 
not  exist  and  the  motion  is  predicted  to  be  an  up  and  down  motion  which 
is  shown  in  Calculation  IV  and  V  in  Table  15.  The  calculated  barriers 
are  illustrated  in  Fig.  3. 
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TABLE  15 
Comparison  of  Barrier  for  N2F2 

No.                      Ref. 

Parameters 

Etot.    (a.u.) 

ZNNF 

Rnn 

RnF 

I 

C2v  (Cis) 

109.98 

1.2197 

1.2U10 

-78.87U0723U 

II 

C2h  (trans) 

10U.09 

1.2218 

1.2UU1 

-78.86889607 

III 

0.5(1  +  II) 

107.035 

1.22075 

1.2U255 

-78.871Ii8U21 

IV 

C2(\ J 

111;.  07 

1.2208 

1.21+26 

-78.7582l;577 

V 

Cgw  ( ) 

180.0 

1.2119 

1.1M, 

-78.5U333813 

„        .         E(II)  -  E(IV)       E(II)  -  E(IV)             E(I)  -  E(IV)       E(l)  -  E(V)  _ 
barrier     =  6o#^6                s  20U.356                    =  72.7057            =  207.657      7mote 

Exp#       E(II)  —  E(V)  =  32  kcal/mole   (37) 
Cal'       E(II)  —  E(V)  =  8U  kcal/mole   (36) 

-78.1* 


C5-78.6 


X 


^-71 


-79.0 


Fig.  3  Barrier  for  N2F2as  function  of      XNF 
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V-7.  Oxygen  difluoride  (F20) 

For  the  F2O  molecule  with  C2V  symmetry  (RqF  =  l.UlA°  (13)),  the 
lowest  energy  configuration  is  calculated  to  be  bent  with  anJj^(Sf  =97.16° 
(at  fixed  R^=1.5>0A°)  as  compared  with  the  experimental  angle  of  103.8° 
(38).  The  bond  length  R£p  (at  angle  of  97.16°)  is  calculated  to  be 
1.1792A0.  Also  the  bending  and  stretching  force  constants  are  calcu- 
lated to  be  1.1U97  x  lO"11  erS/rad?  (observed  "fezFOF  =0.55  x  10"11  erg/rad2) 
and  h0.9k  x  10*  ^/cm  (observed  ^kofi   =5.57  x  io*  dyne/cm  (lit)). 

It  is  of  interest  to  compare  the  orbital  energies  as  a  function 
of  internuclear  angle  with  that  publisted  by  Buenker  and  Peyerimhoff 
(39).  The  calculated  results  are  illustrated  in  Fig.  h   and  these  were 
compared  with  c alculat ions  by  Buenker  and  Peyerimhoff  (using  full 
LCAOSCF  theory)  also  in  Fig.  U.  The  comparison  in  Fig.  Ij.  shows  a 
satisfactory  agreement  except  for  5a^,  Ib-^  and  3bg .  It  is  noted  that, 
in  full  LCAOSCF  calculation,  the  lowest  doubly  degenerate  orbitals 
(due  to  Is  electrons  of  two  Fluorines)  and  the  next  singly  degenerate 
orbital  (due  to  the  Is  electron  of  Oxygen)  were  neglected.  Therefore 
the  comparison  is  valid  starting  with  the  third  level  in  the  Fig  given 
by  Buenker. 
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V-8.  Hydrogen  Peroxide  (H2O2) 

For  the  hydrogen  peroxide  molecule,  the  calculations  were  performed 
to  find  the  barrier  and  the  equilibrium  configurations  with  R^u=0.97A°, 
ZH00=100°,  R~=1.U°A°  and  dihedral  angle  of  93.6°  (13).  Unfortunately, 
the  present  calculation  does  not  show  the  existence  of  a  dihedral 
angle.  The  result  of  the  calculation  shows  that  there  exists  a  Cis 
and  trans  configuration  for  H2O2  just  as  for  ^^2'     Tlie  structural 
parameters  are  illustrated  in  Fig.  £. 


Fig.  5  Structural  Parameters 
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V-9.  Homonuclear  diatomic  molecules 

For  the  diatomic  molecules,  the  equilibrium  configurations  of  H2, 
Li2,  N2  and  F2  were  calculated.  An  attempt  was  made  to  find  a  new  set 
of  optimized  parameters  by  changing  the  Slater  exponents  and  beta 
zeroes  (/J  ).  The  equilibrium  bond  lengths  and  force  constants  are 
listed  in  Table  16.  It  is  noted  that  B2,  Cz  and  O2  molecules  are  excluded 
in  the  calculations  because  they  are  ground  state  triplets,  and  the  Be2 
molecule  ip  non-bonded.  The  necessary  partial  derivatives  were  found  and 
listed  in  Table  16. 


TABLE  16 

Equilibrium  Configurations  and  Some  Partial  derivatives 
for  homo  nuclear  diatomic  molecules 

Re  (A°) 

^rOo5^) 

average  partial  between (S 
and  also  ^3°    ±  1 1n 

.E)o±  0.1 

Cal. 

Obs. 

Cal.    Obs. 

C^sV 

(%Ms 

(*%V 

(^k_ 

H2 

0.7U6 
(0.89) 

0.7U 

10. 1U 
(6.5) 

5.6 

-0.5U 

-0.3h 

19.18 

1.0555 

Li2 

2.179 
(2.78$ 

2.67 

0.81* 

(o.Ui) 

0.25 

-2.9I1 

-0.022 

3.07 

0.52 

Nz 

1.1U 

1.09 

52.16 

23.10 

-0.1*7 

-0.01 

61.06 

2.3U 

F2 

- 

1.12 

l.iiU 

56.1*1 

U.U5 

For  those  molecules,  an  attempt  was  made  to  find  a  new  set  of 
optimized  parameters  fitting  experimental  equilibrium  bond  lengths 
and  force  constants.  The  Slater  exponent  and  yO   were  considered  in 
adjusting  the  parameters  to  optimize  agreement  for  calculated  potential 
surfaces.  The  procedures  used  to  find  a  new  set  of  Slater  exponents 
and  beta  zero  ( $    )  were  as  follows : 
First:  Change  the  Slater  exponent  by  ±  0.1  and  beta  zero  by  ±1.0, 

then  find  the  calculated  equilibrium  bond  length  and  force  constant. 
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Second:     Find  the  partial  derivatives 


"3R 


(filtfl 


•3*r 


"3/5° 


Third: 


where  "S%  stands  for  Slater  exponent. 

Find  the  new  Slater  exponent  and  by  following  equations: 


as  follows 


(S.i) 
(5.2) 


then  approximate  the  dR,  ds  and 

dR  =  *R  =  R«p.-  Roil. 

mS  =  ^S  =  Saew  -  S&'- 

Fourth:     Solve  equation  (5.1)  and   (5,2)  to  find  the  new  Slater  exponent 
and     A      .     For  example: 
(l)  H2  molecule: 


(5.3) 


-0.06 

.-U.83 
therefore 


(2)  Li 


2  case: 


f-OJtfl 

-0.59 
therefore 


(3)  Ng  case: 
-0.05  " 

-29.06 

•  4 

therefore 


'-5.U     -3.U 
19.18    1.055 

Snew  -  0.923 
ytfnew  =^.157 

2.9U       0.022 
3.07       0.52 

Snew  =  0.U9 
/?new  =  8.59 


-0.U7 
61.06 

Snew  =2.87 
fine*  =  -9.5 


-0.277 
0.1^57 


-0.16U 

-o.Uio 


^S 


-0.01 

0.87 

'*5~ 

2.3U 

4 

-35.5 

Aj6° 

(5.h) 


(5.5) 


(5.6) 


hi 


From  these  results,  (5>.h)and  (5.5)  are  consistent  but  ($,6) 
is  completely  inconsistent  because  the  force  constant  is  too 
high. 
Fifth:   Finally,  repeat  the  calculation  using  the  above  new  set  of 
input  parameters  (Snew  and  LJ   new).  The  result  for  the  H2 
and  Li2  molecules  are  listed  in  parentheses  in  Table  16. 

VI.  CONCLUSION 

The  Approximate  Self-Consistent  Field  M.O.  theory  was  applied 
to  the  calculation  of  electronic  structure  for  several  selected  mole- 
cules and  these  results  were  compared  with  results  by  full  LCAOSCFT 
calculations,  if  they  were  available.  The  most  important  features  of 
these  calculations  are  to  predict  equilibrium  configurations  (bond 
length  and  bond  angle)  and  the  barrier  by  the  variation  of  internuclear 
distances  and  angles  using  the  ASCFT.  The  EHT  calculations  are  limited 
tc  hydrocarbons  which  give  good  results  (refer  to  calculation  of  Ethane 
molecule).  In  general,  the  barrier  predicted  by  ASCFT  calculation  is 
in  satisfactory  agreement  with  the  experimental  value. 

An  attempt  was  made  to  find  a  new  set  of  input  parameters  (Slater 
exponent  and  ^S°   )  to  improve  the  prediction  of  potential  surfaces 
for  certain  homonuclear  diatomic  molecules  but  the  attempt  failed 
except  for  the  H2  and  Li^  molecules  (see  Table  16). 

The  features  of  Potential  Surfaces  of  molecules  can  be  predicted 
with  the  approximate  self-consistent  field  theory.  Qualitatively,  the 
agreement  between  calculations  and  experimental  observations  is  quite 
good.  It  is  conceivable  that,  with  the  adjustment  of  certain  parameters 
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such  as  resonance  integrals  (  fl°  ),  Slater  exponents  and  penetration 
integrals,  quantitative  agreement  with  experimental  results  can  be 
obtained . 
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I       Computer  program  for  Approximate  Self-Consistent  Field 
Theory  (ASCFT) 
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APPENDIX 


II   Computer  output  of  Water  (H2O) 
(Sample  calculation) 
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